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ABSTRACT 


For a wide range of block sizes and earthquake magnitude, 
the probability that a free standing rigid block will survive with- 
out overturning when subjected to earthauake ground acceleration is 
determined. Survival statistics are obtained by comouter simulation 
which consisted of calculating the block's motion when it is sub- 
jected to each of a large number of different artificial earthquakes 
and determining for each earthauake whether the block overturns or 
not. Two tvpes of artificial earthquake are used and their prop- 
erties and the method of generation of each are discussed. From 
comparison of the results obtained with each, it is concluded that 
the simpler earthquake mode!, which consists of impulses occurring 
at random time intervals, can safely be used to estimate survival 


probabilities. 
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NOTAT ION 


mechanical energy 

ensemble average of x(t) 

D'Alembert force impulse 

D'Alembert force 

transfer function of a linear differential equation 


= +! when angular displacement is positive 


-| when angular displacement is negative 
height to center of gravity of block 
= one half the height of the block 
moment of inertia 
imap 
i 
subscripts 
kinetic energy 
mass of block 
number of impulses in random impulse model earthquake 
approximate white noise 
potential energy 
scale factor for filtered white noise model earthquakes 
block size parameter 
semi-diagonal of block 


correlation function 
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discrete values of filtered white noise acceleration 
velocity reduction factor due to eneray loss at impact 
power spectral density 

relative velocity response spectrum 
theoretical undamped response spectrum 
spectrum intensity 

time 

time interval 

period of free vibration 

ith impulse in random impulse model earthquake 
impulse magnitude 

weight of block 

semi-width of block 

impulse magnitude 

time average of x(t) 

base acceleration 

relative displacement of oscillator mass 
filtered white noise 

slenderness ratio 

random numbers 

Dirac delta function 

fraction of critical damping 

angular displacement of block 


time interval 
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db = non-dimensional displacement parameter 


i) circular frequency 


Other symbols are defined when used. 
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CHAPTER | 


INTRODUCT ION 
l,l Introduction to the Problem 


In designing structures for earthquake resistance it is 
common practice to apply a static horizontal load in place of the 
actual inertia load. In the case of tall slender or top heavy 
Structures, it is possible that the overturning moment due to the 
eich seismic load is greater than the opposing moment due to the 
weight of the structure necessitating special design of found- 
ations and pilings to prevent overturning. However since the 
actual seismic load is dynamic rather than static and the over- 
turning of such structures seldom occurs during even severe earth- 
quakes, it is of interest to consider the dynamic response of 
these inverted pendulum type structures to earthquake ground 
acceleration. Two examples of this type of behaviour that have 
been reported in the literature are: (a) petroleum cracking 
towers in California during the Arvin-Tehachapi earthquake of 
1952, and (b) golf-ball-on-a-tee water towers during the Chilean 
earthquakes of 1960. These structures survived with relatively 
little damage when other more stable appearing structures were 
severely damaged. It was observed that the cracking towers in 
California had stretched their anchor bolts which suggested that 
their behaviour had been that of a rigid body in a rocking motion 


rather than that of a cantilever beam undergoing elastic vib- 


ration. 
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In 1963 Housner (1)* considered the problem of a rigid 
block resting on a horizontal base which was accelerated by random 
impulses. He found that tall slender blocks when subjected to this 
type of excitation were more stable against overturning than they 
appeared to be when analyzed under the application of a static 
horizontal inertia load. Housner also stiolisd Tay fy WO 
geometrically similar blocks of different size were subjected to 
the same ground motion, the larger block would be more stable than 
The smaller... In considering the probability of,survival of a 
block subjected to the random impulse excitation, Housner was only 
able to determine whether this probability would be greater than 
or less than 50 per,cent. That is, he derived a formula relating 
the size and geometry of the block to the magnitude of the 
excitation to give a 50 per cent survival rate. 

An earlier study of a similar problem was described in 
a paper by Muto, Umemera and Sonobe (2) in 1960. This paper 
describes the dynamic testing of models of slender structures 
resting on a base which was given horizontal sinusoidal motion. 

It was found that the harmonic motion necessary to overturn 
structures would have to be of very long periods and of large 
amplitudes. Such motions do not occur in earthquakes and it 
was concluded that tall slender structures would not be as 


easily overturned by earthquakes as might be supposed. 
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* Numbers in parentheses refer to references in Bibliography 
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To determine the probability of a rocking structure over- 
turning under earthquake action is an example of the first passage 
time problem of random vibration theory. There exists no exact 
theoretical solution to determine the probability that a random 
process will exceed specified limits for the first time ina 
given time interval. Several approximate solutions have been 
derived but to apply them to a particular process it is necessary 
to know certain statistical propertiesof that process (14,15,15). 


\f the random process that is being considered is the response 
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of a linear time invariant system, it is possible to determine 
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the required response statistics 
and the impulse response function of the differential equation 
by applying superposition techniques via the Duhamel integral. 

It will be seen that the differential equation of 
motion for the case of a rocking block is not continuous but that 
it is step-wise linear so that superposition cannot be used. 
This makes it impossible to calculate The statistical properties 
of the angular displacement of the block from those of the 
excitation using standard analytical techniques. In such cases 
response statistics can be estimated from data obtained from 
repeated observations of the system. With the speed of modern 
computers, it is possible To obtain a large number of observations 
by simulating the system on a digital computer. This is the 
procedure that is used in the investigation described in this 
thesis. The details of a block's motion are computed to 
determine whether or not it overturns when it is subjected in turn 


to each of a large number of different earthquakes of the same 
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magnitude. From the number of earthquakes that the block survives 
without overturning, the probability of that block Surviving an 
earthquake of that magnitude can be determined. 

There does not exist a large number of records of 
actual earthquake ground motion of a given magnitude that could 
be used as the excitation in such a procedure. It is necessary 
to generate random functions suitable for use as models of earth- 
quake ground motion. The two types of earthquake models used in 
this investigation were: (a) ground acceleration consisting of 
impulses which are random in time, and (b) ground acceleration 
defined by a segmentally linear function consisting of randomly 
varying discrete values at small time intervals At connected by 


straight lines. 
LZ Aimsiiof the Thes&s 


|. To confirm and extend the results obtained by G. 
Housner (1) concerning the probability of survival of a rocking 
block subjected to earthquake acceleration modeled as a series 
of impulses arriving at random time intervals. 

2. To obtain similar results using a more comprehensive 
model of earthquake acceleration obtained as a form of filtered 
white noise and compare with those obtained above. 

3. To use simulation techniques, modeling the system 
on the 1.B.M. 360/67 digital computer at the University of 


Alberta Computing Center, to determine the required statistics. 
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CHAPTER || 


ARTIFICIAL EARTHQUAKES 
Z.1 Random Impulse Model Earthquakes 


Ze clvat Reenonee Spectrum and Spectrum Intensity 
A concept which has been developed to investigate 
and delineate the character of earthquake ground motion and to aid in 
determining structural response to earthauakes is the response 
spectrum. The response spectrum of an earthauake is the curve 
obtained when, for a series of linear oscillators whose undamped 
natural periods of vibration range from .10 seconds to 3.00 
seconds, the maximum response of each oscillator is plotted 
against its period for a particular damping factor. The response 
function which is commonly used is the relative velocity. Consider 
the spring-mass-damper system shown in Figure 2.1] which is excited 
by the base acceleration x(t) and exoeriences the displacement 
y(t) of the mass relative to the base. The equations determining 
the relative velocity response spectrum value, Sy(n,w), are: 
y + 2nwy + wy = -x(t) 

Sy(n,w) = |y|max. Ze 
in which n is the fraction of critical damping and w denotes the 
natura! circular freauency. Repeated calculations of S, for 
different natural frequencies and different damping factors for 
a given acceleration record X(t) produce a set of curves such 


as are shown in Figure 2.2 for one component of the Taft, Calif. 


earthquake of 1952, (|). 
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Figure 2.3 shows the response spectra curves obtained when 
the spectrum values from eight actual earthquakes are averaaed. 
Though the response spectrum can be used to determine the expected 
maximum response of linear structures to a given excitation, for 
purposes of the present investigation, average spectra will be 
used as the means to compare artificially generated earthquakes 
Tesdetermine their.suitabibity. 

Because the structures in an area which is experiencing 
an earthquake will have various natural frequencies, it is 
desirable to have some average over all vertinent frequencies 
which would be an indication of the severity of the earthauake. 

A single parameter which accomplishes this is the Spectrum 
Intensity, Sly, which is defined as the area under the velocity 
response spectrum with damping n between the periods .10 
seconds and 2.50 seconds. 

To form the average velocity spectra of the actual 
earthquakes (Figure 2.3), the eight individual spectra were 
first normalized so that the undamped spectrum intensities were 


all the same and then averaged. 


2.1.2 Random Impulse Model of Earthquake 
One of the first attempts to explain earth- 
quake accelerations in terms of a random process was made by 
Housner in 1947 (7) when he showed the similarity between re- 
sponse spectra for actual earthquakes and response spectra for 


accelerations consisting of impulses arriving randomly in time. 
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Figure 2.1 Damped Linear Oscillator Subjected To 


Base Acceleration. 
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Figure 2.2 Typical Velocity Response Spectra. 


brs) ALO 7. Si 
Natural Period - Seconds 


Fiqure 2.3 Average Velocity Response Spectra For Recorded 
Strong Motion Earthquakes. 
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Consider again the linear system shown in Figure 2.1 and 

the corresponding differential equation: 
y + 2nwy + RAY = -x(+) 

Let the base acceleration consist of an impulse of magnitude X 
oecurring at itime t = 0, x(t) = X6(+), where X is a real constant 
and 6(t) is the Dirac delta function. 

Considering the base, such an acceleration impulse is 
equivalent to a step change in the base velocity. This may be 


shown as follows: 


x + Ax te 
Ax = f Cie th oC cy gee. 
x 26 
Considering the mass, m, the base acceleration may be 

replaced by a D'Alembert force: 

f(t) = -mx(t) = -mX6(t) 
From impulse -momentum considerations 

f(t)dt =m dy 


The mass m thus experiences a step change in its relative 


velocity due to the impulse in the base acceleration as may 


be shown: 
€ +e 


+ 
=f.) £Ct) dt eae -x6(i dt = 
y =E m se 
Thus an acceleration impulse of magnitude X which causes 


a step change in the velocity of the base is experienced by the 


mass as a step change in velocity of magnitude X. If the system 
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were at rest prior to the impulse, the effect of the impulse 
would be the same as free vibration of the mass with initial 
conditions y(0) = 0 and y(0) = -X for which the solution for 
y(+), the relative displacement ,is 
y(t) = Esme sin wt 
For a base acceleration consisting of a sequence of 
N imoulses of magnitude v; occurring at random times tT; where 


heel ,2,.+..N, the resulting motion could be found by super- 


position to be: 


ne e NWt-Ti) sin w(t-t,) Zee 
| W 


Model earthquakes composed of random ground 
acceleration impulses or, what is equivalent, random step 
changes in ground velocity were produced on the computer by 
using a subroutine for generating pseudo-random numbers 
(17). The duration of the model earthquakes was chosen to 
be 30 seconds which has been estimated to be the approximate 
duration of the strong-motion phase of a very large earth- 
quake (10). It was assumed that 100 impulses occurring at 
random time intervals for a total elapsed time of approx- 
imately 30 seconds would be .a suitable model to simulate the 
earthquake ground accelerations. Thus the mean time interval 
between impulses is required to be 0.30 seconds. This was 
obtained by selecting the time interval as a random number 
from a uniform distribution between 0 and 0.60. The impulses 


were taken to be of constant magnitude, Av, but random in 
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sign. To determine the sign (direction) of the impulse, a 
random number was selected from a uniform distribution between 
-|.0 and +1.0 and the sign of the impulse set the same as the 
sign of the random number. In this way there is equal 
probability of any given impulse being positive or negative. 
Thus 100 times, t;, separated by intervals of random length 
and 100 impulses, v; = £ Av, were selected to define the 


model earthquake. 


2.1.5 Theoretical Response Spectrum 
To determine the theoretical undamped 
velocity response spectrum value, denoted) S$. 4(0,0);, for’ this 
excitation according to Equation 2.1, we must consider the 


maximum value to |y| for n=0. From Equation 2.2 with n=0 


yet t Wen 1) Tete 
j= 
Differentiation yields 
N 
y = - Vi Os wT = 1) 
i=| 
N 
= ~-. V; (cos wt cos wt; + sin wt sin wt;) 
i=| 
N 
SVC. iv (COS wt) cos wt 


N 
ote Vie oikn wt ;) sin wt 


The addition of these two harmonic functions results ina 


single function with amplitude: 
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Substitution into Equation 2.1! and expanding yields for the 
theoretical undamped spectrum value corresponding to natural 


frequency w: 


| ee 


N > N N 
OU a) et ot 2 Vivy Cos EL Goreme 1d) 
i=! i=l jaitl : 
But vane + Av so the first summation becomes 
N 
Fl ye SEN Cm 


In the double summation, the product vV;v; has equal probability 
of being either +(Av)2 or -(Av)%. Due to the randomness of the 
times t; the cosine will have equal probability of being any 
value between -] and +I so that each Term in the double 
summation has equal probability of being any number between 
+(Av)4 and -~(Av)?, With the large number of terms involved 
SOG ei eras +2 +1 for N = 100), the double summation would 
be equal to zero and the spectrum value becomes: 


(0,0) =|NCavy2 


=[N av 


= 10 Av Za 


VF 


for the case of [00 impulses. 


Since the spectrum value, is independent of w, the 
undamped response spectrum curve is a straight line. When 
this model of earthauake acceleration was first considered (7) 


the average zero damped spectrum curves of actual earth- 


iVuh ae \ 
PO UMAY te Taal ee 


‘Ae foe ie Wy, ke ee ISAS Wie i rit (6 A ata ALR A 4 ha 
i , by K ren " ise ~ ; : si iin Th 


= 

ae 
ae 

—v 


aia heii dee tibia << 
| rt Cae a 
a0, isbn WA zy a & i 
' a) } t. 
a | 
GAB cece Pee by yiert ice % p teitiae ern wie 
| ia 


Keay der! it gi BA yi ft heres? ao 


a a 
‘: — g . : : . Pink : 
veo no ™“paye fri ait ee ee 
f PT ave 
' 
= j 
( 
4 = 
Jib VA i 1” 
ry pa} 
PT TAL ale se iil t ik ey at 
A : , 7 ¥ 
Mt Lene ne) % it a | fete Ww) ee 
emis We WET Ty Fata) en ilia 
mee ety ee r see, ean a TM Ee 
' 4 ’ 
on) ' 


Te Ls i + ojrhey, j re PER a . au ifabsé, | ore es jit 
‘ bf i” am i 
CAN ete) Mo oo Bete eee of ie Ly - tne? 


Hila ey ied UT oid othe «iy ld 49% pa TR wd ee i 


ol 4 ) 
a} pean 4 au taey hu rere HOR i eee: or tad 
—- rie. \y fi : 
: oe ; 7 F . 
' m vane : 4 af % 
i Le : 4 4 Ae , ie fl 
vi AS 
j : ‘ A o ( 
» Ram why * 
j J th vas wi : : y ani a 
eae | Li | i ? | apenas 
e po . ri ain . a . 


a 
: 7 


quakes that were available for comparison had been calculated by 
a mechanical analogue that incorporated approximately one per 
cent of critical damping. This was sufficient to eliminate 
the hump from the zero damped spectrum so that the Straight line 
spectrum of the random impulse type acceleration compared 
favourably with the actual earthquake results and it was 
assumed that this type of ground acceleration could approximate 
earthquake action. 

For vibration periods in excess of |.5 seconds, the 
average undamped spectrum as shown in Figure 2.3 becomes a 
flat curve independent of period so that for structures with 
large natural periods, the random impulse tvpe of excitation 


could be a good model for earthquake around motion, 


2.1.4 Computer Evaluated Response Spectrum 

To determine the actual velocity response 
spectrum curves for a random impulse model earthquake gen- 
erated as outlined above, the linear system of Figure 2.| was 
simulated on the computer. The simulation consisted of 
calculating the details of the motion of the oscillator mass 
due to the base acceleration x(t) defined by the earthquake's 
impulses v; at times tj. By repeating the simulation for 
oscillators with different natural periods and different 
damping and determining the maximum relative velocity attained 


in each case, the necessary data to plot the spectrum curves 
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was determined. 

In section 2,1.2 it was shown that the only effect of 
a ground impulse was to cause a step change in the velocity of 
the oscillator so that between impulses it moves as free vib- 
ration, Thus the motion of the oscillator was calculated as a 
series of free vibration problems defined as: 


yj(t) + 2nwy(t) + wy, (4) = 0 


fOr (OP Sin KS At, 

with ¥iisG0 hes y, pdt, 1) 

and y; (0) = y,-, 4t._)? + Vv; 
MEeneeTiie SUDSSHILD hy hee ely Z piaien-a 100 denotes the time interval 


in which the motion is being calculated (i.e. the interval 
between ur and te and At, = Tay aitpitstimhe lengeth of this 
time interval. 


The motion of the oscillator is started by the first 
impulse, v;, of the earthquake so that the initial conditions 
for the first interval were yi, £0? = 0 and y (0) = Wine 

The displacement and velocity were determined at The 
end of each interval to be used for the initial conditions of 
the next interval. The velocity was also calculated at all 
times when the value of the displacement was zero. 

Whenever there was a zero crossing within an interval, 


a new initial condition problem would be considered in the sub- 
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interval from the time at which y = 0 to the end of the interval. 
The initial conditions for this sub-interval would be y(O) = 0 
and y(0) = velocity at the time of the zero crossing. 

For small damping the general solution for the dis- 
placement in all cases was taken as: 

y(t) = Y e7RWT sin (wt + A) 
with the constants Y and \ determined from the initial cond- 
itions. 

From all the values of velocity calculated during the 
oscillator's motion, the maximum velocity is determined and 
this is the velocity response spectrum value for the part-- 
icular value of natural period and damping factor being con- 
Sidered. The response spectrum curve is formed by plotting 
the spectrum values for natural: periods of .10, .20, .30.., 
2.50 seconds and joining them by straight lines. Because of 
the random nature of the excitation the response spectrum 
for any one earthquake is not as significant as that obtained 
as the average for a large number of earthquakes. Average 
spectrum curves were obtained by averaging in each case the 
results obtained in analyzing 150 random impulse model 
earthquakes with a given impulse magnitude. 

The average response spectra obtained for model 
earthquakes with impulse magnitude Av = 0.25 and damping 
ratios of 0, 0.02, 0.05, and 0.10 are shown in Figure 2.4. 
Typical undamped average spectrum curves are shown in Figures 


2,5, 2.6, and 2,7. Also shown in these figures are the 
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average ordinates of the calculated spectrum curve and the 


gerresponding straight line spectrum, Sts predicted by theor- 
etical considerations. Table 2.1! summarizes the results ob- 
tained in the calculation of average undamped velocity spectra 
of random impulse model earthquakes with impulse magnitudes 
Banoo trom 0,125 f1./see.. to 0.55, fte/sec. 

The spectrum intensities were calculated as the 
area under the spectrum curve and the average ordinate de- 
termined as that value which a straight line spectrum must 
have tn order to have the same spectrum intensity. 

In comparing the calculated and theoretical undamned 
spectrums, it was found that the calculated spectrum was 
independent of the period of vibration as predicted, but that 
the average ordinate of the calculated spectrum was con- 
sistently higher than the predicted spectrum value, 

S4(0,w) S/N Av 
for all impulse magnitudes considered. When corresponding 
values of the average ordinate of the calculated spectrum 
and the theoretical spectrum value were plotted as shown in 
Figure 2.8., it was found that there was a linear relation- 
ship between the two characterized by the slope: 


average ordinate = 
9 
vt 


A corrected expression for Sy in terms of the 
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number of impulses and the impulse magnitude which would 
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agree with the average ordinate of the calculated spectrum would 


be 
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Table 2.1 


Comparison of Theoretical and Computer Evaluated Spectrum 


Values For Random Impulse Model Earthquakes. 


| | i i 
Impulse ; Theoretical | Calculated Spectrum | 

Magnitude Spectrum value | i | 

| S4(0,0)=/N - Av | Average Spectrum | 

| Ordinate intensity 4 

| Bate le 

| |.25 | 3.36 

1.50 | 4.03 

i 

: 1.87 5.04 

2.00 | 5.39 

2.50 | 6.78 

50 | 9.50 
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Figures 2.5 - 2.7 Average Undamped Velocity Response Spectrum 
Curves For Random Impulse Model Earthquakes. 
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Figure 2.8 Theoretical Undamped Response Spectrum Value Versus 


the Average Ordinate of the Calculated Spectrum for 


Random Impulse Model Earthquakes. 
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2.2 Filtered White Noise Model Earthquakes 


2.2.1 Definitions and Concepts From Random Process Theory 
A random process x(t) is one which cannot be 

described by a single record, for as the process is repeated, 
the records obtained continually differ from each other. It is 
necessary to consider the entire ensemble of possible time 
histories x, (1), Xo(t), Seana x (1) from which a statistical 
description of the process can be derived. The first order 
probability distribution describes the distribution of values 
x(t,) which occur in the ensemble for t=t, by means of a graph 
of the first order probability density function p(x, t)). 
The fraction of ensemble members for which x(t) ) lies between x 
and x + dx, which is equivalent to the probability of finding 
x(t) between x and x + dx in any given record of the random 
process, is given by p, (x,t) dx. The expected or mean value 
of x(t) is given by E[x(t))] = Hs xp (x, 7) ) dx where 
E{ ] indicates an average taken across the ensemble. Higher 
order probability distributions describe the distribution of 
of products x(t x(t), x(t) )x(t5)x(t2), etc. across the 
ensemble and provide a measure of the mutual relation between 
the values of the random function at different times. 
Sufficient probability information for most processes can be 
determined from the first and second order probability dis- 


tributions. The second order probability density function 
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is denoted Po(xy ty; Xot,) and Palxy,tlsx2,t2)dx) dx. is the 
probability of finding any one of the sample functions 
x, (fF) in the interval (x) ,x)+dx)) at time leat and in the 


interva| (x., X> + dx.) at “hime. f.= to. 


A stationary random process is one for which the 
Statistical parameters do not vary with time. For such a 
process the first order probability density function is time 


invariant, bp (x,t) becomes Bp) (x) and the exoected value of 


| 
XPamany) Timers *alx}%= je. xp, Oddx. The interdependence of 
of the values of x at two instants of time depends only on 

the time interval between them, Tt = t, - Ty, so that the 
second order probability density function is a function of 

TH Polxy, Xo t) Because the parameters do not vary with time 
it becomes reasonable to consider averages taken with respect 
to time. According to the ergodic assumption for a station- 
ary random process, statistical properties determined by 
averaging over time for any one sample will be equal to those 
determined by averaging across the ensemble. For a given 


sample function x(t) the mean value of x as determined by 


the time average, denoted ees)? 


1 
eetiam: ly fo exe 
* Witoass2 fT 
By the ergodic assumption every sample function will have the 
same average which will be equal to the average across the 


ensemble: x = E[x] 
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The variance, o , is the ensemble average of the 


square of the deviation from the mean: 


gS ELK — E[x])29] 48f oieGes Epar pl Geax 


| 
For processes with mean zero, i.e. E[x]=0, the variance is 
identical with the mean square: ae = E[x*]. 

A statistical measure involving second order prob- 
ability densities is the correlation function defined as the 


average of the product x(t)x(t+1) and denoted RMT) 


co co 


Bete = AARC Seg if f X)X5P(X) »Xo,7) dx, dx. 


In terms of the time averages the correlation function is 
defined as 
ny uF 
Retry = cert eye iii [Pe Oe ex CR ede 
Tee EL! @ ok 


A frequency analysis of a stationary random process 


can be accomplished by considering the power spectral density, 


S(w), which is defined in terms of the correlation function 


by using the Fourier integral transform: 
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For a stationary random process S(w)dw is a 


measure of the average power available in the frequency 


ee 


interval w tow + dw. Whereas the response spectrum discussed 


earlier (Section 2.1!) is a concept used only in earthauake 
engineering, power spectral density is used more generally 


as a means to describe the frequency content of a stationary 
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random process. A random process for which the power spectral 
density has a constant value over the entire range of frequencies, 
i.e. S(w) constant for all w, is referred to as white noise. 

A normal or Gaussian random process is one for which 
the first order probability density function is a normal curve 
and all higher order probability densities are similarly 
functions of known form which can be defined in terms of the 
power spectral density. In this paper we wil! be concerned 
with Gaussian random processes with zero means and the re- 
sponse of a linear system to this type of random excitation. 

When a linear time invariant system is excited by a 
zero mean Gaussian random process, the response will also be a 
zero mean Gaussian process. The power spectral densities of the 
excitation and response functions are denoted Sia (w) and 


Soyt(w) respectively and are related by 


eC S, tw) rg 


Sout (w= 
where G(iw) is the complex frequency response or transfer 
function of the linear system. CG(Ciw) is a function of fre- 
quency such that when the excitation input to the system is 
elt the response is bite ue 

2.2.2. Filtered White Noise Earthquake Similation 
The relation between the excitation and 
response spectral densities (Equation 2.4) forms the basis for 


the method used by Housner and Jennings (10) to produce 


Gaussian stationary random functions with a specified spectral 
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density suitable for use as models of earthquake ground accel - 
eration. If a linear spring-mass-damper system is excited by 
white noise so that Sin (w) = constant, then the magnitude and 
form of the spectral density of the response will depend on the 
transfer function of the system. By proper selection of the 
resonance and damping characteristics, the system can be de- 
signed to produce a random motion with the desired statistical! 
properties. The output obtained from such a system is ref- 
erred to as filtered white noise. 

Using a linear system such as ithat shown in Figure 2.1 
as the filtering system, with the base acceleration, Gon being 
the white noise excitation and the absolute acceleration of the 
mass, z(t), being the filtered white noise response, Housner 
and Jenninas produced an ensemble of eight pseudo-earthquake 
accelerograms which they proposed to use in studies of the 
earthquake response of structures. The same approach was used 
in the present investigation to produce an ensemble of fifty 
mode! earthquakes. 

For the filtered white noise to be used as earth- 
quake acceleration, it must have the same spectral density 
as the recorded ground accelerations of actual earthquakes. 
However, the statistical properties of earthquakes are known in 
terms of the average velocity response spectra shown in Figure 
2.3. In their paper, Housner and Jennings show that the 


relation between the velocity response spectra, S,(n,w,t), 
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and the power spectral density, S(w), of recorded earthquakes 


is’given by: 


Sy (nit) = 1.796 /xStwd , (1 - eW 2Net/ 2-44, Ds 
2nw 


where t is the time duration of the earthquake ensemble. For 
zero damping this reduces to 
$(0,w,t) = 1.150 frStwt 
When the above relation was used to determine the 
spectral density to correspond to the average response spectra, 
it was found that a simple algebraic function of a form which 


is convenient for linearly filtered white noise could be used 


Oto 'C0)) 2 
0201238 (1 + Sa) 
S(w) = aE BL ZG 


a 
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2.2.3. Generation of Filtered White Noise Model 


Earthquakes 


Consider again a linear mass-spring-damper 
system as shown in Figure 2.1. The equation of motion can be 
written as: 

y + 2bcy + cMy = -x(t) 
where y is the displacement of the mass relative to the base, 
x is the forced acceleration of the base and b and c are system 
constants. If z is defined as the absolute displacement of the 


mass so that z = y + x then this equation may be written as: 
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To determine the transfer function, Gliw), between 
x and z ( which is also the transfer function between x and 3 ) 
we set x = e!4T ang z = Gliw)elt jp Equation. 2.6. Solving 


for Gliw) yields: 
em 2bw 
€ 


2 
| ~ oS 4 ; aby 


C4 C 


Glin) = 


lf the base acceleration x(t) is white noise its 
power spectral density will be constant, say S; fw) = a and 
The absolute acceleration of the mass, z(t), will be filtered 
white noise with a spectral density Pe eae determined by 


Equation 2.3 to be: 


1 
is 
-—_- 
E 
a) 


Sout fo) = 


If the filtered white noise obtained from such a 
system is to be used as a model for earthquake ground accel - 
eration, it must have a spectral density as given by Equation 
2,0. Joes ing Sout Se? from Equation 2.8 equal to S(w) from 
Equation 2.6 the required values of the coefficients of the 
differential eauation of the filtering system and the magni- 


tude of the input spectral density were found to be: 


c2 = 242, b% = 0.410, a = 0.01238 = S, Cw) 


ng Bhi A 
nate nee Si 
A: : L. 


7. af Ap. : 
7 ! 


ash. Wee WENN A 
ie oh) rae OO q 
oty | oun! i 
hy i 4 yy 
0) 4 : Lat : 
"4 hone 


a) ; 


i, aos vee My her 
HI i abut, PaaS ment te att aes ae 


12 bie e fabwt ed ry it oiut SeNeRIe et: ate a sti 9 


Dae a 


ph vs Oe Waa MLNS A nila “ea sti ie 


; og | i 
4 ' | fe 
oe, 4 : al niin i a 
ood ‘ x i ih ; 
oy Pepiple- nae - ty Hi Tal : ' 
ee ae | 


att. waled wii 40 Cecilie pees Seen eee 4th. ae 
, A , ¢ -4 ’ ' Hl 
a ld bys eed 4 bigtenoe ea) 4 tnd yh a ‘ey es : 
, ’ , s sy ' 
Sie) a ie 2 Ok i iy he ha ervleedae 
‘ mY, , ey ie 
4 


MO DeMHATOW (Ch). peer daned Jet page, A nw salon 
7 a 


en une 


cen i Pi . pe mitt tips) vide m aks the nts bi aa a By, 7 
+, ey i) . o. Ops ol ra an ee , 7 Na 
eee Land 9) eat oh, aheotatnen “ot siete ; r 7 be ft 

a i ue, et? ye ue fh, 


were 


aa Pi: ys 


ey 


27 


To use this system to produce model earthquakes, 
white noise was required for use as the excitation. An 
approximate form of Gaussian white noise, denoted by N(t), 
was formed by generating a sequence of Gaussian random 
numbers (17) which were used as the ordinates of the function 
separated by time intervals At and connected by straight 
lines to form a segmentally linear function. It was shown by 
Housner and Jennings (10) that if At is taken as 0.025 sec. 
and the random numbers have a mean of zero and a variance 
9% = eh ., then the function N(t) will have a spectral density 
close to unity, Sy fw) = ).0.+ @05. for thesranae of, freq= 
uencies of concern. However, it was required that the input 
acceleration x (+) have ayspectral density oe Gu e= aca AS 
the spectral density is a second order function of the process, 
Mieracceleration was taken as x(7) =ja N(t) to yield the 
spectral density of x (+) as 


ie me 2 
S 5 6w) = fa ) oR a? a 


as required. 


Consider again the differential equation of the 

filtering system in terms of the relative displacement y, 
y + 2bey + c“y = -x = —/faN(t) 2.9 

To determine the time record of the absolute acceleration 


of the mass, Z (+), .which forthe vadues of the coefficients 
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determined above will be the required filtered white noise 
accelerogram for a model earthquake, this equation must be 
integrated to evaluate 

2(t) = y +x = -2bcy - oy. 
To start the integration initial values of the relative 
displacement y(0) and velocity y(0) are required. Because the 
excitation was assumed to be a stationary random process, 
y(t) would also be stationary and it was sufficient to choose 
the initial values from the probability distributions of 


y(t) and y(t). It can be shown that their variances are 
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lf 8 and y are two numbers taken at random from a Gaussian 


distribution with mean zero and variance unity then 
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will ensure that these initial values are members of Gaussian 


processes with mean zero and the required variances. 


The integration of Equation 2,9 was accomplished on 
the digital computer using a fourth order Runge-Kutta method. 
The value of > (+) was determined at 0.025 second intervals 


and these values joined by straight lines to produce the model 


7 i ae am PAS) ae 
, : rr 


a * 
ey 4 M 


en 
hs 
ny ae ie 
| mH : ieee bier 


af 


eh Se dk 
ball 


ive) Qe. Weis 


Ne heeded i eae Mies i hal oer ~anoss Mallee. wed 
} , te ; 


easy aaeome ey Be a bins: cab ca. Pe + one 


oe 


an ste bean ibe Pit ety . aA: At doitirignyn me: 2 


haw wey Bs 
ae ww Hartt a baths 
M in ey “i Bayt ay any i pert 
ery ne ‘ x mY “i, whats si ie 


ak =) = { 


es) 


earthquake accelerogram. To obtain accelerograms of thirty 
seconds duration, 1200 random numbers were required to produce 
the white noise excitation and Similarly 1200 numbers were 
determined to define z(t), 

By repeating the integration, each time with a 
different sample of the white noise excitation, as many 
different accelerograms as required can be produced. However, 
because the computer time required was considerable, only 
fifty acceleroqrams were produced to serve as the excitation 
ensemble in the present investigation. 

The filtered white noise accelerograms were modified 
slightly to enhance their suitability as realistic models 
for earthquake ground acceleration by ensuring that the 
ground velocity tends to zero at the end of the earthquake. 

To achieve this, a corrected base line defined by the function 
Co + Cy ateet os was superimposed on the accelerogram and the 
velocity of the ground found by integrating the acceleration 
as measured from this corrected base line starting with zero 
initial conditions on the velocity and displacement. The 
constants Co» Cy, and C. were determined so as to minimize 

the mean square of the velocity and then the acceleration 
ordinates were determined with respect to the corrected base 
line to form the model earthquake accelerograms. The 


integration to determine the ground velocity is an easy and 
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exact calculation because of the acceleroaram being seamental ly 
linear. 

2.2.4 Response Spectra of Filtered White Noise 

The relative velocity response spectra of 
the model earthquakes were determined by considering again a 
linear oscillator whose base is excited by the adjusted 
accelerograms which will be denoted z(t). The equation 
determining the spectrum value for a particular freauency w 
and damning ratio n are: 
y + 2nwy + oy = -z(t) 


S(n, w= |¥ | max 


The integration of this equation was carried out 
using a Runge-Kutta method starting with zero initial 
conditions on the relative displacement and velocity. The 
integration interval used was half the basic interval of the 
excitation so that the relative velocity of the oscillator 
was determined at .0!125 second intervals and the calculated 
velocity with the largest magnitude was taken as the response 
spectrum value. For each model earthquake spectrum values 
corresponding to oscillator periods of 0.10, 0.20, ....1.00, 
h.20, 1.40 1...2.80 and 3.00 seconds and damping ratios of 
0, 0.02, 0.05 and 0.10 were determined to completely specify 
response spectra and the spectrum intensities. From the 


spectrum values for each of the fifty mode! earthquakes, 
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average spectrum values for the ensemble were determined to 
define the average response spectrum curves. 

To judge the suitability of the model! earthquakes, 
if is preferable to compare their average response spectra 
with average spectra of recorded earthquakes. However, the 
model earthquake spectra cannot be compared directly with the 
recorded earthquake spectra given in Fiqure 2.3 because they 
are not of The same duration and in Equation 2.5 it is shown 
that the spectrum values depend on the duration. The average 
duration for the recorded earthquakes was 18.6 seconds whereas 
the duration for the pseudo-earthquakes was 30 seconds. 
Equation 2.5 can be used to determine the expected average 
response spectra for the mode! earthquakes. 

Figure 2.9 shows plots of the response spectra 
formula and the average response spectra of the model earth- 
quakes. The agreement between the two is better than that 
obtained by Housner and Jennings with their ensemble of 
eight model earthquakes as would be expected due to the 
larger ensemble. Such agreement speaks well for the 
consistency of this method of producing pseudo-earthquakes 
as well as showing that the accelerograms produced should 
prove satisfactory as mode! earthquakes. 

The value of the zero damped spectrum intensity, 

SI was determined for each of the fifty model earthquakes. 
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The values obtained ranged from Slo = 2.68. To Sl = 95,04. 
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The average for the ensemble was ST = 3.13 and this value was 
taken as the parameter to indicate the strenath of the ensemble. 
To use the ensemble to represent earthquakes of various 
strengths, the acceleration values defining the fifty earth— 
quakes were multiplied by a scale factor, say 9, to produce 


an ensemble with average spectrum intensity ou ee ae e ts 
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CHAPTERD | 
ROCKING BLOCK EXCITED BY RANDOM IMPULSE MODEL EARTHQUAKES 


3.1 Free Vibration of a Rocking Block 
Consider a rigid block as shown in Figure 3.1 of 


weight, W, mass, M, moment of inertia, I about an axis 


0? 
through the corner 0, and with the angle a between the 
diagonal and the side of the block defining its pronortions 
and referred to as the slenderness ratio. When such a 
block is given an angular displacement, @, measured from the 
| vertical as shown, then providing 8 < a, the weight will 
cause a restoring moment about 0 of magnitude WR sinla - 6) 
where R is the radial distance along the diagonal from 0 
to the center of gravity. Taking clockwise as the positive 
direction for angular displacement, the equation of motion 
for positive displacement can be written as 
oF 6 2 -WResin(a@.- 6), 02 6 2a 

For tall, slender blocks such that the angle a is small 
enough to use the approximation sitnto? =16)5S 0 =09,% this 
equation becomes 

at 8 ~ WRe = -WRa Sih 
A property of the block, denoted a, which is analogous To 


the natural circular freauency of a harmonic oscillator mav 


by defined such that q2 = WR For a uniform solid block 


Io 


Io = 4 ur? = 4 WRA so that a? = 54 , Because the value 
aoe. q AR 
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of q depends only on the lenath R it will be referred to as the 
Size parameter of the block. Dividing Equation 3.| Dy fy, 
yields 6 - 976 =)=q5a. 
In terms of a non-dimensional quantity 9,defined 
as ¢ = , the equation of motion becomes 
6 - 9% = -q? OS gs] ae 
For *initial sconditions ¢ t= ég and 6 ie ire i iy eae 0 
which correspond to the block being given an initial 
displacement 6, = 694% and released, the solution of this 
equation may be shown to be 
otach sullequgey Boshyag hy) eOlSioesal Bn3 
The above derivation is applicable only for ® 2 0. 
Thus Equation 3.3 describes the block's motion only as it 
rotates back into the vertical position. The time t, at 
which the block would be in the vertical position can be 
determined by substituting ¢ = 0 in Equation 3.3 to yield 


il = | 
ty = £ wCosh {-—— ) O08! 


q -%, 


If there is no energy loss during the impact at 


6 0, the block will rotate about the corner at 0! until 
@ = -0, and then fall back to the vertical position. To 
describe this phase of the block's motion, it is necessary 


to determine a different equation of motion and its 


solution subject to initial conditions of zero disnlacement 


and known initial angular velocity at t = t;. The equation 
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of motion for negative displacement can be written as 
I,6 = WRsinla+6), -aSoa0 


With q and @ defined as above, this equation becomes 


o - a*4 = a2, <li Spd oe 375 


Foneinitraltconditions = 0 and b b6 art = t,, the 


solution of This equation ts 
% 6 
Ops olb =. (Sinh- grat) 4 Coshmatt-F1), 3.6 
: mire: Ge ke 
Equation 3.6 would describe the block's motion as it 
rotates from 6 = eee and back to the vertical. 

Each subsequent half-cycle of the block's motion 
could in turn be described by either Equation 3.2 or 
Equation 3.5 depending on whether the block was moving on 
the positive or the negative side of the vertical. 

For the case of negligible energy loss, the block's 
motion would be symmetrical about the vertical position and 
the time t+, as determined by Equation 3.4 would be equal to 
one quarter of the period of free vibration, i), so that 

Tees Conn eset} 


| 
I-o,, 
Thus the time period of the vibration is dependent 
on the amplitude. 
For the case when the energy loss at impact is not 
negligible, the angular velocity of the block after impact 


would be less than that just prior to impact and a velocity 
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reduction factor,r 5 1.0, could be defined such that 


are  TWwheRS 6 is the anaqular velocity before 
6| ?| 


impact and 8. is the angular velocity after impact. Due to the 
loss of energy, the amplitude and hence the period of each 
half-cycle would be less than those of the preceding hal f- 
cycle. The free vibretion of the block then would consist of 
initial oscillations of large amplitude and long periods 
with amplitude and period decreasing for each subsequent 
half-cycle until the block comes to rest in the vertical 
position. 
3.2 Overturning By Random Impulses as Per Housner 

As it was shown above, the motion of a rocking 
block is stepwise linear with the discontinuity between 
steps caused by the necessity to alternate between two diff- 
erent equations of motion. In considering random vibration 
problems, the statistical properties of the response 
process can usually be determined from those of the excitation 
process by means of a relationship determined from the 
differential equation describing the motion. For the case 
of a rocking block, there is no single differential 
equation of motion and such a relationship cannot be deter- 
mined which would apply continuously so that it is not 
possible to apply the usual methods of statistical dynamics. 


Thus a statistical description of the block's motion cannot 
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be determined from the properties of the random impulse ex- 
€itat ion: 

Housner (1) in considering the problem of a rock- 
ing block subjected to random impulses was able to produce 
some measure of the probability of the block overturning by 
considering the average energy input due to the impulses 
in an analysis similar to the following. 

Consider a block such as is shown in Pagure ls. 2 
with the base on which it is standing being subjected to a 
random impulse model earthquake excitation, x, consisting of 
impulses, v; = + Av, occurring at random time intervals. The 
effect of a horizontal acceleration impulse, v;, would be 
equivalent to an impulsive inertial force, F,At, acting through 
Miececenter Of Gnaviryeot the (ohock such ahak 

Bek Mw a ea AN, en Paro RA Ae 
where the subscript ij denotes the effect of the ith impulse 
and M is the mass of the block. This force impulse. causes 
a moment impulse about the corner 0 of magnitude 
F,At R cos(a-6). Limiting consideration to tall slender 
blocks for which a is small enough that cos(a - 6) may be 
taken equal to unity and from impulse momentum considerations 
of the rotational motion of the block about the axis 0 we 


obtain 
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Thus the step change in the angular velocity of the block due to 
the ith ground impulse may be determined as: 


AG. 
Ty te 


Sy 


The kinetic energy of the block would undergo an increment change 
due to the ith impulse of 
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The average change in energy per impulse is found by averaging over 


all impulses to obtain 


(v.)4 


But there is equal probability of any one impulse being positive or 


negative so that the average over all impulses would give ve = 0 
2 


and since for any impulse (v,) = (AYE, then Cae = (Av) 


Thus the average increment in energy becomes 


If there were no energy loss during rocking, the energy of the block 


would increase so that on the average, the block's energy at the end 


of N impulses would be 


(MR)2 
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As the enerqy of the block increases so does its amplitude of vib- 


ration. The maximum amplitude the block can attain without overturn- 
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ing is a and to vibrate with this amplitude the energy of ane Dlock 
must be equal to the change in potential energy as the block moves 
pom thcvertical position to 6 = a, that ts 

P= WRC] - cos a) = 3 WR a4 3.10 
From Equations 3.9 and 3.10 it can be determined that on the average, 
the condition for the block to be vibrating with sufficient eneray 
to attain an amplitude a after N impulses is 
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When Equation 3.11 is satisfied, then the probability of 
the block overturning would be equal to the probability of the 
next impulse being in the direction to increase the block's eneray. 
As there is equal probability of each impulse being in either dir- 
ection, Equation 3.1! represents the condition for 50 per cent 
probability that the block will overturn. 

[In chapter two it was shown that for a random impulse 


model earthquake the theoretical undamped response spectrum 


would be 


S4(0,w) = aval 


In terms of the response spectrum, Equation 3.11 may be written 
as 
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ape MR2i4 5 


For a particular block subjected to a around motion 
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with a given response spectrum, it is only possible to determine 
from the above theory whether the probability of overturning is 
greater than or less than 50 per cent. 


For a uniform solid block Lae : MR? and Equation 3.12 


S460,0) = af #.Rg : 3.13 


relating the size and proportions of a block to the response 


becomes 


spectrum of the earthquake for which it would have a 50 per 
cent probability of overturning. This relationship shows that 
for constant a the larger the value of R the stronger the 
earthquake that may be survived. That is, for two blocks of 
the same proportions, the taller block would be more stable 
than the shorter. This can also be explained by substituting 


into» Equation: 5:/ttozobrain 
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which shows that the angular impulse AA, imparted to a block 


by a ground impulse v; is inversely proportional to the size of 


the block. 


3.3 Survival Probabilities Determined by Simulation 

To test Housner's theoretical results and to develop 
moré complete data on Phe sprobabl Way of a block remaining up- 
right when subjected to earthquake ground motion modeled as 
random impulses, the system shown in Figure 3.2 was simulated 


on the digital computer for various values of the slenderness 
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Fatio, a, the block size parameter, q, and the earthquake mag- 
nitude parameter, Sys Each individual simulation consisted of 
calculating the details of a given block's motion when subjected 
to a particular random impulse model earthquake of a given mag- 
nitude to determine whether or not the block would overturn. 
From the number of survivals out of 1000 such simulations, each 
with a different earthquake, the probability of the block not 
overturning, termed the survival rate, was determined for that 
particular block and earthquake magnitude. By varying the para- 
meters, results were obtained for the significant range of 
block sizes and earthquake magnitudes. 

In each simulation, the excitation, x, consisted of 
a different random impulse model earthquake with a duration of 
30 seconds generated as outlined in Section 2.1.2. The un- 
damped velocity response spectrum for these earthquakes was 
foundwro%be S,(0,0) = : . 10.Av so that earthquakes of varying 
spectrum value were produced by varying Av, the magnitude of 


the impulses. As shown in Equation 3.7 an acceleration imnoulse 


Av causes a step change in the anaular velocity of the block 
Z 
5 = BM oe 


: WR = 39g os 7 
where as previously defined q? ies ie . In terms of the non 


dimensional displacement function ¢ 6/a , the step change in 


velocity is ; 
ee 3.14 
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Equation 3.14 shows that blocks which have the same parameter q 
(i.e. the same semi-diagonal R) but have different slenderness 
ratios a will experience the same magnitude of angular impulse 
Ad as long as they are subjected to model earthauakes with im- 
pulse magnitude Av such that Av/a is the same for all blocks. 
Since the block's motion in terms of $ depends directly on the 
magnitude of the impulse wen it was to be expected that the 
survival rate of a block with given value of q would depend on 
the ratio Av/a or since the response spectrum is defined in 
terms of Av alone, on the ratio Sy/a. Thus by determining the 
survival rate for a particular ratio S,/a and various values of 
the size parameter q, the coordinates were obtained to plot a 
curve of survival rate versus q. A number of such curves for 
different values of S\/a were obtained by means of which the 
expected survival rate for any particular block and earthquake 
were shown. 

In each simulation to determine the motion of a 
given block due to a particular earthquake, the total time of 
the motion was divided into a number of intervals within each 
of which the block's motion was determined as in initial cond- 
ition problem. The effect of each ground impulse of magnitude 
Av was to cause a step change in the angular velocity, or an 
angular impulse of magnitude Ad. Between impulses the block 


moves as in free vibration so that its motion consists of a 


consecutive series of free vibration problems, each one 
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starting just after the occurrence of one impulse and ending just 
before the next impulse. The initial conditions in each case 
would be the final conditions of the previous interval plus the 
effect of the impulse. That is, the initial displacement would 
be the displacement at the time of the impulse and the initial 
velocity would be the velocity at the time of the impulse plus 
the step change in velocity due to the impulse. 

In considering the free vibration of a rocking block, 
it was found that two equations of motion were required, one 
which held for positive displacements and one for negative dis- 
placements of the block. This meant that in addition to a new 
initial condition problem being initiated by each impulse, it was 
also necessary to consider a new problem each time the block 
passed through the vertical (zero displacement) position. The 
initial conditions for this case would be zero displacement and 
initial velocity equal to the velocity just prior to impact 
multiplied by the velocity reduction factor, :r,ito eccount for 
energy loss at impact. If energy losses were considered neg- 
ligible r could be taken as unity. 

The two equations of motion 
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by defining H as H = -! when 6 2 0 


H | when 630 
so that H would change value each time the block passed through 
the zero displacement position. 


The initial condition problems were then defined as: 


6,(t) - q* o;(t) = Hq? = for OS + < At; 3.15 
and 6,(0) = $;_,(At;_)) + ad; 
or 6; (0) = d; Catz er 


where the first velocity condition applied for an interval in- 
itiated by an impulse, in which bd; was the step change in vel- 
ocity due to the impulse and the second condition for an interval 
initiated by passage through zero displacement, in which r is 
thes velocity reduction,;factor.+ Ihexsubscriph.~iadenotes the in- 
terval and At; is the length of the ith time interval. 

In the generation of the model earthquake to be- used tn 
the simulation, the first ground impulse, v;, was always taken as 
being positive, v, = tAv, and subsequent impulses randomly selected 


as positive or negative, v; = + AV, eas deseribeduine PeCt WonuZ al. 2 


Each ground impulse will cause a step change in velocity, Ad;, where 
. q? Vj 
ip = So fener 
g a 


The first ground impulse, v,, will start the block moving with an 
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initial angular velocity Ady. setting v,; positive ensured that in 


the first interval the block would experience positive displace- 


ment so that H = -I| and the initial conditions for the first inter- 
val were 
e a? 
$60) = 0, 610) 2 grees eer 
g a 


In each interval the solution for the displacement was 
taken in the form 6 = | + AGt Te notte meunere she constants A 
and B were determined from the initial conditions. From the 
ratio A/B it was determined whether there was a possibility of 
a peak value, ¢ = Onan, OG a 2ero Crossing, 0 =)0,) before “tne 
next impulse. If no zero crossing occurred, then the interval 
would be terminated by the next impulse and the values of 6 
and $ at the time of the impulse were calculated and used to 
_ determine initial conditions for the next interval. 

As the block would overturn if the anale it makes 
with the vertical becomes greater in magnitude than a 
because the moment due to the weight would then be acting to 
increase the angle rather than restore the block to the 
vertical position, the necessary conditions for survival is that 
lo|<1.0 throughout the block's motion. Whenever the displace- 
ment was calculated, either at the end of an interval or as the 
peak value in an interval, a check was made to see if Tits cond= 


ition was satisfied. 


The calculations of the details of the block's motion 


iy i 


iit ; 


Runa 7 wae ro 
ba hee SBNT: bernie ay LE Lda | i ei 
cw : ' ~~ 


Moa hea i Awa Pease yp a ty tae Si on wer Ving nee 


r ci fi } tery 5 end «oe ay | 14 Aon ab f ‘wt une i ee | Tt icles 


i volt 4 BS het (ate iu 
r ‘ 


Se tiovns Hyon’ Ihr bir ym g eC eh bible wilt 1 


f 

i ‘ "ao 

TO To OAM HAIRS 4 hy Sata ete med 3 4 
Ls ee AAs 

/ j iT 


TO. hy een “$41 Sadibert uo OSES mee it ~~ | 
wir rs C1) 220% ieee ie Cguwe® aa | ie i 
ie) 


f ne} aye w 2 ss ai, te - ; 
aw viet 1h ae “ir Tl t ‘cans aves: cap. a4 oR ui 


Auibibie an ert iri4) ad yt Sather thts AIK Mi pe i al é 


- ‘ A 
ee Ser ee ee ae aht tobrbtsdal 008 aif 


Oe oy oh MR hse igo un ina | 
Pas het ob ph pe stir. aK rere 
Ue ae ANTNG 38 ee ais aw sah ia ene a 
7 nN ie tf ie CY SAN ent sy saan sets ion eae eae 
ie oat ene PP) mbhvine tor conrad vores! of sot nm i 


' oy ati oer’ les a i bas Bidsald aby. ‘tom 


47 


were continued until a value of displacement was determined such 
that [o|71.0 or until the last impulse of the model earthquake 
had been accounted for with either a failure or survival being 


recorded for that particular earthquake. 


3.4 Results of Computer Simulation 
In the curves of Figure 3.3 the survival rate is plott- 
ed against the block size parameter q for various values of the 
ratio S/o. These curves present the data as determined by 
computer simulation for the ideal case of perfectly elastic im- 
pact for which the velocity reduction factor, r, is equal to one. 
For each particular value of S\,/a the curve extends 
from a minimum value of q such that for q 5 q,;, the survival 
rate will be 100% to a maximum value of a such that for aq : 
Gee tne survival rate will be zero. Between these two ex- 


tremes as q increases the survival rate decreases. For a un- 


4. so that increasing gq 


\ 


corresponds to decreasing R. A curve of R verses q is plott- 


iform block it has been shown that q 


ed in the top part of Figure 3.3 and superimposed on it are 
some sketches showing the sizes of blocks with a = 0.10 
corresponding to particular values of q. It is apparent that 
these curves show agreement with the general result predicted 
from consideration of Equation 3.13 derived by Housner's 


approximate theory that a taller block would be more stable 


than a shorter one of the same provortions. This can be dem- 
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onstrated by considering two different blocks, one 100 feet high 
by 10 feet wide and the other 25 feet by 2.5 feet, each subjected 
to an earthauake WITH, = 5.5 47. sec. The slenderness *ratiosin 
each case is a = 0.10 so the particular survival curve of inter- 
est is S\/a = 35. As shown by the dashed lines in Figure 3.3 the 
probability of the taller block surviving would be 79% compared to 
a probability of only 11% for the shorter block. 
Equation 3.13 predicts theoretically the response 


spectrum value for which a block would have a 50% survival rate. 


ee 
In terms of the size parameter, q ae , this equation becomes 
g 
Se pee 3.16 
a q 


The graph of Equation 3.16 is shown in Figure 3.4 where S,/a is 
plotted against q. Corresponding to the various values of S\/a 
in Figure 3.3, the value of q was determined for which the sur- 
vival rate was 50% as determined by the computer simulations. 
These pairs of values were plotted as additional points in 
Figure 3.4 to show the agreement between the theory and the 
simulation results for this particular case of 50% survival 
probability. 

Thus the limited results predicted by the approximate 
energy approach for the case of a random impulse excitation 
were confirmed by the data obtained from the computer simul- 
ation of the same situation. In addition, the simulation 


analysis gives results for a much broader range of circum- 
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stances. 

Each survival rate value used in the data of Figure 3.3 
was determined by subjecting a particular block to each of a 
randomly generated set of one thousand model! earthquakes and de- 
termining how many earthquakes out of the sample set it would 
survive without overturning. If a different sample of one 
thousand earthquakes had been used, it would be expected that a 
different survival rate would result due to the statistical 
nature of the problem and the limitations imposed by a finite 
sample size. 

To obtain an indication of the statistical consistency 
or "goodness" of the simulation results, one particular case, 
S\,/a = 25, was repeated seven times, each time with a different 
set of model earthquakes. The variation in the survival rates ob- 
tained is shown in Figure 3.5 where the range of values obtained 
for certain block sizes is shown superimnosed on the original 
curve for this case, S,/a = 25, from figure 3.3. From being neg- 
ligible near the extremes of zero and 100% survival rate, the 
variation increases to a maximum in the area of 50% survival rate. 
fomades 1.50 the ayes rates obtained ranged from 44.3% to 
52.1% with a mean value of 47.7%. There was good general agree- 
ment among the results obtained from the eight different sets of 
earthauakes with no outstanding inconsistencies. A larger sam- 


ple size would yield more reliable survival rates but to use 
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more than one thousand earthquakes would cause the required computer 
time to be prohibitive. 

To illustrate the effect of eneray loss due to impact, 
the survival data for S\/a = 25 was determined with the velocity 
reduction factor having the values r = 0.95 and r = 0.90. The 
curves thus obtained are also plotted in Figure 3.5 for com- 
parison with the curve for the ideal case, r = 1.0, from Figure 3.3 
to demonstrate that the effect of any eneray loss is to increase 
the probability of a block's survival. 

lt would be difficult to assess what value r might have 
for any actual situation as it would depend on the material prop- 
erties of the ground and the structure. For the case of complete- 
ly inelastic impact application of the conservation of angular 
momentum yields an expression for r in terms of a, the slender- 
ness ratio. Consider Figure 3.6 which shows a block just prior 
to impact rotating about the corner 0 with angular velocity 8 
and angular momentum about 0! given by (Io - MR=)6 | + MRO | 
(h cos a - w sina). After impact the block will continue its 
rotation about the corner 0' with velocity 85 and the angular mo- 


mentum I, . 65. Equating angular momentum about GF iusto or 


to impact to that just after impact yields 
156) - MRO,[R - (h cosa -wsina)]=Ij162 3.17 


But R-(h cos a - Ww Sin a) = 2wsin a = 2R epics and I,:= Ip 
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from the geometry of the block so Equation 3.17 becomes 


Io6) - 2MR* sin2a8, = Ty85 3.18 
For slender blocks such that sin asa the velocity re- 


duction factor is obtained from Equation 3.18 to be 


iG SS 6/82 | 
4 
subst turing 1-3 Z MR? fer a uni tormeapliock yields sihe 


relation 


showing that the energy loss varies directly with the square of 


jhe stenderness ratio. 
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Figure 3.2 Rigid Block Subjected to Base Acceleration. 
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CHAPTER IV 


ROCKING BLOCK EXCITED BY FILTERED WHITE NOISE MODEL EARTHQUAKE 


4.1 Introduction 

To obtain a different appraisal of the probability of the 
rocking block surviving without overturning in an earthquake, the 
system of Figure 3.2 was simulated on the computer using filtered 
white noise model earthquakes for the ground acceleration x. 

The ensemble of fifty standard earthquake accelerograms 
produced as outlined in chapter two was used for the excitation. 
Each earthquake was defined by 1200 values of the acceleration at 
intervals of .025 seconds, connected by straight lines to produce 
a segmentally linear function. To vary the severity of the earth- 
quake ensemble, the acceleration values were multiplied by a scale 
factor Q. 

For .a given strength of earthquake, i.e. a given value of 
Q, a block with a particular slenderness ratio a@ and size para- 
meter q was subjected to each of the fifty earthquakes in turn. 

For each earthquake the motion of the block was calculated to de- 
termine whether the block would overturn or not. From the number 
of survivals out of the fifty earthquakes, the survival rate was 


determined for that particular block and earthquake magnitude. 


4.2 Forced Vibration of a Rocking Block. 


Consider a rigid block with properties a, R and W as 
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defined in chapter three and shown in Figure’ 3) .freé“tovrockon 
the ground which is experiencing acceleration x(t) as shown in 
Figure 3.2. The effect of the ground acceleration can be account- 
ed for by an inertia force, se x, acting through the center of mass 
of the block as shown. Taking clockwise as the positive direction 
for Bea tiay displacement, 6, and x positive in the direction shown, 
the equation of motion for positive displacement is 


1,0= -WR-sinla - @) + : x R cosa - 0).) O8S0 “Sa 4.) 


For tall slender blocks, the approximations sinla - 6)= a - 86 


and cos(a - 8)= | may be used and with q@ = re as defined in 


chapter three, Equation 4.1 becomes 


Neh iain Alera: 


In terms of the non-dimensional auantity @ defined as 


¢ = = the equation of motion for positive displacements becomes 
Pebcepeku lee, 142 1308 Gre (phen 4.2 
eas aS : 


Because of the change in the direction of the restoring 
moment due to the weight, a different equation describes the motion 
when the block is experiencing negative displacement. This 
equation is: 

16 = WR sinta FO) + - xR cos(a + 6) Sy 465.0 
In terms of q and $ this equation becomes 
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lf the: block, initially at rest, is’subjected:to a positive 
base acceleration, x > QO, it would tend to experience a positive 
displacement and Equation 4.2 would apply. When $ = 0,) Miecan ibe 
seen from Equation 4.2 that a necessary condition to obtain a 
positive angular acceleration (4 > 0) and hence start the block ro- 
tating is x > ga. That is, the base acceleration must be greater 
than the product of the gravitational acceleration and the slender- 
ness ratio. If the ground acceleration is such that 0 < x < ga 
the block will experience a translational motion with the ground 
but it will not rotate. 

Similarly for the case of negative ground acceleration, 
%160, consideration of Equation 4.3 with ¢ = 0 reveals that for 
the block to acquire a negative angular acceleration, 6 << OE, 
is required that x < =G0 

Thus the condition on the magnitude of ground accel- 
eGearion Peroni to initiate rotational motion of the block can 
be stated as |x| > ga. 

By using the variable H as defined in chapter three, 


the two Equations 4.2 and 4.3 were replaced by one equation 


$ - q%6 = Hq as q? a > eal Scere 4.4 
ga 
with H-= -l. when ¢-2°0 and H= +i whenig-< 0 . 


The solution to Equation 4.4 for the nondimensional 
displacement function $(t) depends on the type of function de- 


fining x(+), the ground acceleration. 
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The particular case of interest is ground acceleration 
defined as a linear function, 
x(t) = a + bt 4.5 
Substitution into Equation 4.4 yields 
b - 976 = Hq? + q@ (atbt) 4.6 
DefinindQse:H a. andy] Se-ecieige 
gga 


ga 


g2(C + Degt) Duy, 


" 


¢ 976 
for which the general solution can be taken to be 
o(t) = Aedt + Be~9T - © - Deqt 4.8 
The coefficients A and B may be determined by the in- 


itial values of the displacement, $(0), and velocity, $(0). 


4.3 Computer Simulation of a Rocking Block Subjected to Filtered 
White Noise. 

The time record of the filtered white noise ground 
acceleration consists of consecutive time intervals of length 
At = 0.025 seconds within each of which the acceleration is a 
linear function. Within the ith time interval, (i = |,2,..1200) 
which extends from tT; = Gt HW Araror tee iAt the acceler- 
ation was defined by values R; and Reay at times tT; and Ti+ 


respectively and expressed as 
(Ri+1 TRG) 


x; (tf) = Rj - At 
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with t as the time measured from the beginning of the interval so 
“HAG sli) oy Oe aaa go 

The motion of the block was calculated as a series of 
consecutive initial condition problems with the displacement and 
velocity at the end of one interval serving as the initial cond- 
itions for the next interval. Using Equation 4.7 with C and D 


défined from comparison of Equations 4.9 and 4.5 as C = H + = 


and D = Ss , the initial condition problems were defined 
Atqga 
as: 
ce Z a aa < < 
Pell iG a CD) SG tCor OgT) FORO TS AT 4.10 
and 6; (0) = $;_) (At) 


As was shown in Section 4.2 no rocking of the block will 
occur until the magnitude of the ground acceleration becomes 
large enough. If the initial value of the acceleration satisfies 
the condition IR, |>ge, the proper value was assigned to H 
rier LG ui) 20, ee lala R, <0) and tne motion tn the Tigst.in= 


terval would be determined from Equation 4.10 with starting 
conditions $(0) = $(0) = 0. 

| f IR, |<ga , then it became necessary to determine the 
first interval within which |x| becomes equal to ga. If R 5+ 


was the first value of R; for which Ra [>ga, then motion of 
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the block would begin in the jth interval at the time Tog. ber srnet 
interval when [xj (tg ]=ga. Measuring time, t, from the start of 
motion, then within the subinterval of length, At' = At - t, from 
the start of motion to the end of the jth interval, the accel- 


eration would be x (4) = a + (Rut i-Ro) Piha 
nae 


in which R, is the vaiue of ground acceleration when motion 


starts. If the ground acceleration is positive, Ro. woe. and 


fee liand {he constant Cin Equanwtonea. 7. would ie 
ea re hae rl Be 
ge ga 
lf the ground acceleration is negative, Baie econ ee rim] 
3 ; ee (Rj4}-Ro) ; : 
enamagatn,.U. = 0. Veniring 243s | ne mon (Onn ihe 


At qga 
first sub-interval was determined from 


! 


6 C1) - g°6 j(t) = q“Dat , Che cher 
with $500) = $)(0) = 0. 


In each interval or subinterval, the solution for the 
displacement was taken in the form 
6(t) = Aedt + Be~St - C - Dat 4.11 
with the constants C and D as defined in terms of the excitat- 
ion and the constants A and 8 determined from the initial cond- 
itions. Differentiation yields the corresponding expression 


for angular velocity: 
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lf the displacement attained a value such that ale if.) 
overturning of the block would occur. To check for failure it was 
necessary to know the magnitude of the displacement at any maximum 


or minimum peaks that occurred in an interval. 


Peak values of ¢ will occur at times t, when $(t,) = 0. 


Reference to Equation 4.12 shows that +, can be determined from 
qty ae aoe qt) oN 
- 762 - D. Multiplying through by e and defining a 
fies , 2 
new variable S =e yields the quadratic equation 0 = As - Ds - B 


the roots of which are given by S = = (Di wRDF + akeeoe Ux 


0 = Ae 


\ 


(p* + 4AB)<O no peaks occur in the interval. If (pé + 4AB)>0 the 


displacement function given by Equation 4.11 has two peaks. 

lf it was determined that one or both of the peaks occurred 
at times t,<At, then the peak value(s) of @ were calculated to 
determine if overturning would occurr in the interval. If there 
were no peaks in the interval, the value of ¢ at the end of the 
interval, o¢( At ), was calculated. 

Whenever the value of ¢ was calculated at a peak or at 
the end of an interval, it was necessary to check that the algebraic 
sign of ¢ had not changed compared to the previous value of $ 
calculated. A sign change would indicate that a zero value of $¢ 
had occurred in the interval. Whenever ¢ passes Through zero 
(which corresponds to the block passing through the vertical 
position) the sign of the parameter H changes so that the ex- 


pression for o(t), Equation 4.11, is only valid up to the time 
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to at which o(t5) =. 1) 

To solve for the time t> when o(t,) = 0,.an iterative pro- 
cess, the Newton-Raphson method, was used to determine the root oe 
the equation, 

Oa=sqet amet per Thise nore Dat, 

Uf o, were the last value of the displacement calculated 
before the zero and $5 the first value after the zero so that $ | 
and o> were opposite in algebraic sign, a first approximation to 
to was made by determining the time To at which a linear function 


from $ to bo would be zero. Subsequent approximation to t5 were 


calculated by the sequence 
o(t,) 


Totnt Py Bsr = 
¢ at ates 


with the iteration being continued until |6(t5)[<107° then this 

value of T, was taken as the time at which the zero occurred. 
Following a zero eeoee ing at Time to in the jth Time: in- 

terval, a new initial condition problem was defined for motion 

in the sub-interval from t5 to the end of the jth interval. The 

value of H would be changed to the negative of its previous value. 

Measuring time, t, from the zero crossing, then within the sub- 


fatervaleof length At’ = AT = to, the acceleration would be 


y (Ree eRe 
x(t) =R,+—HL 2 - t 
J At! 


in which ip is the value of around acceleration when > = 0. 
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The constants C and D were defined as 


Cita A) a: Ro 

ga 
(R.,,-R,) 
and Die Seamer 
At'aga 


and the motion in the sub-interval was determined from 
6) = q¢ (4) = q“(C Vik Git) form iOusiehaeinae! 
with o | (0) = 0 
and $)(0) = $(to) +r 
in which $(t5) is the velocity just prior to impact and r is the 
velocity reduction factor as defined in chapter three. 
Theoretically, as stated above, no matter how small the 
angular velocity of the block as it falls to the vertical (zero 
displacement) position, it would continue to rotate after impact. 
This leads to the possibility of the block rocking back and forth 
with very smal! amplitude and smal! periods of vibration, such 
that for practical purposes it could be considered to be standing 
still. To avoid calculating the details of such motion, which 
would involve extensive computer time due to the large numbers of 
peaks and zero crossinas that would occur, a limiting value was 


placed on the velocity before impact and if the velocity was less 


than or equal to this limit, it was assumed that the block came to 


rest. This limiting velocity was taken to be 0.0001 so that if 


é|< 0,000! the block was considered to stop. To initiate motion 
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of the block again would required the ground acceleration to be 
such that Ix |>go so that the process outlined previously would be 
repeated for values of R; in subsequent intervals. 

The calculations of the details of the block's motion 
were continued until a value of displacement was determined such 
that 1 [21 or until the last interval of the pseudo-earthquake 
acceleration had been accounted for with either a failure or sur- 


vival being recorded for that particular earthquake. 


4.4 Results of Comnuter Simulation 

In Figure 4.1 the survival data are presented for blocks 
subjected to the filtered white noise type of model earthquakes 
considering perfectly elastic impact. For various values of the 
ratio Q/a, where Q is the excitation magnitude narameter and a is 
the slenderness ratio, curves of survival rate versus size parameter 
qa, are plotted. 

From comparison with Figure 3.3 i¥ is seen that there is 
good qualitative agreement between the data obtained using the two 
different models of earthquake ground acceleration. 

To make a meaningful quantitative comparison if is nec- 
essary to determine a criterion which relates the theoretical re- 
Sponse spectrum value, S\4, which defines the magnitude of the 
random impulse excitation, to the scale factor Q used to define 


the magnitude of the filtered white noise excitation. 
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Comparison could be made on the basis of the undamped 
spectrum intensity Sl,, the area under the zero damped spectrum 
Curve. “AS Shown in Section 3.1, the natural freauency of a rock- 
ing block depends on its amplitude so during its motion a block 
will have a varying natural frequency. Since Sle provides a mea- 
sure of the earthquake's effect aver a range of frequencies, it 
seems reasonable to use it as a basis for comparison. For the 
random impulse model earthquakes, the spectrum intensity is de- 
termined from the spectrum value as 
Soiree ashe aaa ee, 

For the filtered white noise model earthquakes, the undamped 

spectrum intensity for the standard ensemble is 3.13 so for 

various values of the scale factor Q, the spectrum intensity is 
Steir eat NC 

The relationship between Q and S\4 such that both types of earth- 

quake models would have the same spectrum intensity is 


SAOQU= 2.74 


from which @) O.675 Syt Aatss 
Another possible method of comparison could involve the 
response spectra. Before a block will overturn it must be rock- 
ing with large amplitude and corresponding large natural period. 
For large natural periods the zero damped spectrum curve for the 


standard ensemble of filtered white noise model earthquakes 


levels off at a value of 1.25 so that for a given value of Q 
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Sy = 1.25 + Q. For the random impulse earthquake, the spectrum has 
the same value at all periods, 7 Susicye . EQuating these (values (sc 
that both earthquake models have the same spectrum value corresp- 


onding to long time period yields 


8 


by25>Q = j S yt 


or Q = 0.915 J+ 

It is seen that these two methods of camparison yield 
criteria which agree within five per cent. Using Equation 4.13, 
Goprespondyng values ofthe ratios*s,+/a and 0/a) Tobe used for 
comparisons were determined and are shown in Table 4.1. The sec- 
ond row shows the value of Q/a that should be used according to 
Equation 4.12 and the third row shows the nearest value of Q/a 


available from the data. 


TABLE 4.1 


CORRESPONDING VALUES OF S\+/a and Q/a 


To show the comparison of the corresponding cases of 
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Table 4.1, three pairs of curves are shown in Figure. 4.2, ° In 

most cases, the results obtained using the random imoulse model 
earthquakes are more conservative than those using the filtered 
white noise model. That is, for a given value of the size para- 
meter, q, the probability of survival as read from the Sy+/ curve 
is less than that given by the Q/a curve. As the ratios O/oa and 
Sy+/a decrease, this difference between corresponding cases be- 
comes greater. 

One factor that could cause this difference is the 
possibility of a block not being set into rocking motion by the 
finite accelerations of the filtered white noise. As discussed 
Mmrcection 4.2, the block will not rock untaligthe guound ecce! = 
eration, x(t), attains a value such that Ix|>qa. The ground 
acceleration magnitude depends directly on the scale factor, Q, 
so that decreasing 9 enhances the possibility of the block not 
rocking for at least part of the earthquake's duration. In the 
case of the random impulse excitation, the acceleration attains 
infinite magnitude during each impulse so that the block must be 
rocking during the entire earthquake. Due to this major differ- 


ence, it is to be expected that the impulses are a more severe 


form of excitation. 
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CHAPTER V 


CONCLUSIONS and RECOMMENDATIONS 


5.1 Summary 

Using computer simulation, which consisted of calculating 
the details of the motion of the system on the digital computer, 
the probability of a free standing rigid block surviving an earth- 
quake type disturbance without overturning was determined. Two 
different typves of pseudo-earthauakes, (2a) random impulses and 
(b) filtered white noise, produced as detailed in chapter two, were 
used for the ground acceleration. The data obtained for the case 
of random impulses was presented graphically in Figure 3.3 and 
found to be in agreement with results predicted by Housner's 
analysis (!). The data obtained for the case of filtered white 
noise was presented in Figure 4.1. Comparison of Figures 3.3 and 
4.1! showed good agreement between the two types of excitation 
with regard to the general relationships between survival rate, 
block size and geometry, and excitation magnitude. With regard to 
the actual survival rates, it was found that when comparison was 
made on the basis of equal spectrum intensity that the random 


impulse model earthquake yielded more conservative values for 


the probability of survival. 
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5¢2-Conélusions 

|. Computer simulation provides a useful tool for de- 
termining response statistics for random processes, particularly 
for situations in which the statistical quantities required can 
not be calculated. 

2. Data obtained by simulation confirmed two of 
Housner's results (1!) for the case of a block subjected to ran- 
dom impulses: (a) the relationship between Sy a fand Igri ones0? 
probability of “surviving (Equatiion S96), “and Cb) that TPfPh*two 
blocks of the same proportion experience the same ground impulses 
the larger block will have a higher probability of surviving. 

+. The. probability of survival’ of @ rocking structure 
can safely and conveniently be predicted by subjecting it to 
random impulse ground acceleration with a spectrum intensity 
eaual to that of the expected earthquake. The impulse model 
earthquake is recommended for the following reasons: 

(a) Although the filtered white noise type 

of model earthquake appears to be more realistic, it requires 
considerably more computer time, both to produce the model 


earthquake accelerogram and to calculate the blocks motion 


when subjected to it. 


(b) The ease with which a random impulse 
mode! earthquake can be produced and the relative simplicity 


of calculating the resulting motion of the block makes it 
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possible to consider a larger sample and thus obtain more reliable 
statistical results. 

(c) For the same spectrum intensity value, the 
random impulse model earthquake yields more conservative predictions 


for the survival probability. 


5.3 Recommendations 

|. Most of the survival data was obtained for the ideal 
fee: no enerqy loss at impact between the block and the aground. 
It would be of interest to investigate the type of energy loss 
which may actually occur and its effect on survival. 

2. The model earthquakes used in this paper were of 30 
seconds duration. As the majority of earthquakes will cause strong 
ground motion for less than 30 seconds, the effect on the survival 


rate of varying the duration should be investigated. 
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